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Abstract: By introducing Crossing functions and hyper-parameters I show that 
the Bayesian interpretation of the Crossing Statistics [1] can be used trivially for 
the purpose of model selection among cosmological models. In this approach to 
falsify a cosmological model there is no need to compare it with other models or 
assume any particular form of parametrization for the cosmological quantities like 
luminosity distance, Hubble parameter or equation of state of dark energy. Instead, 
hyper-parameters of Crossing functions perform as discriminators between correct 
and wrong models. Using this approach one can falsify any assumed cosmological 
model without putting priors on the underlying actual model of the universe and its 
parameters, hence the issue of dark energy parametrization is resolved. It will be 
also shown that the sensitivity of the method to the intrinsic dispersion of the data is 
small that is another important characteristic of the method in testing cosmological 
models dealing with data with high uncertainties. 
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1. Introduction 

One of the main goals of physical cosmology is to reconstruct the expansion history 
of the universe and finding the actual model of dark energy. There have been many 
approaches introduced in the cosmological community to reconstruct the expansion 
history of the universe and the properties of dark energy. One can categorize them 
in general to parametric and non-parametric methods [2]. Non-parametric methods 
are least bias by initial assumptions but in many cases they are not efficient and 
estimation of the errors can be also problematic [3-16]. Contrary to non-parametric 
methods, parametric methods are easy to use and they may seem so efficient in 
putting tight constraints on the cosmological quantities but if the assumed param- 
eterization being wrong, which can be easily the case since the actual model of the 
universe is not yet known, then they can be completely misleading and result to 
wrong conclusions. See [17-21] for details of data analysis and methods of paramet- 
ric reconstruction of the properties of dark energy using supernovae data. 

In this paper I introduce Crossing functions and hyper-parameters, motivated 
by the idea of the Crossing Statistic [1], as a tool to discriminate between correct 
and wrong models. It is known that if an assumed model is indeed the correct one to 
describe Gaussian distributed data, then the histogram of the normalized residuals 
should also have a Gaussian distribution, with zero mean and a standard deviation 
of 1. If the histogram instead exhibits significant deviation from the a Gaussian 
distribution, however, then this can be used to rule out the assumed model. The 
Crossing Statistic [1] pushes this well known idea from statistical analysis a step 
further by pointing to the fact that not only should the whole sample of residuals 
have a Gaussian distribution around the mean, but so should any continuous sub- 
sample. In other words we have shown that catching the systematic correlations in 
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the residuals can be used to falsify assumed models. In the current paper I intro- 
duce the Bayesian interpretation of the Crossing Statistic looking for the systematic 
correlations in the residuals by comparing an assumed model with its own possible 
variations. In this approach observable quantities of any assumed model is multiplied 
by a Crossing function and then is fitted to the data and the resultant values of the 
hyper-parameters of the Crossing function decide about consistency of the model 
and the data. In fact derived values of the Crossing hyper-parameters indicate that 
how smooth deviations from the assumed model can improve the fit to the data. If 
improvement in the fit to the data being significant by smooth deviation from the as- 
sumed model, it hints towards the inconsistency of the model and the data. All these 
can be done in the framework of \ 2 statistics to define the confidence limits. Using 
supernovae type la simulated data I will show how efficiently Crossing Statistics can 
discriminate between cosmological models while by using the standard statistics we 
need to assume a parametric form for the cosmological quantities. The core idea 
behind this work is to suggest that instead of parametrizing cosmological quantities 
like luminosity distance, Hubble parameter or equation of state of dark energy, we 
can simply falsify any assumed cosmological model by parameterizing the possible 
deviations from its actual predictions for the cosmological observables. 

In the following I will introduce the Crossing functions and their hyper-parameters 
and argue how they can be used as discriminators between cosmological models. I 
will then use simulated data to show how effectively Crossing Statistics can distin- 
guish among cosmological models. Minimal sensitivity of the method to the intrinsic 
dispersion of the data will be discussed later before I summarize the results at the 
end. 

2. Method and Analysis 

If the real Universe differs from the assumed theoretical model, one would hope that 
it would be possible to develop a statistical test that would be able to pick up on 
this. To these ends we consider the 'crossings' between the predictions of a given 
model, and the real Universe from which the data has been derived. The idea of 
Crossing Statistics is based on the fact that any wrong assumed model and the data 
cross each other at one or two or maybe even more points within the range of data. 
Here I argue that since within the frame work of FLRW universe the cosmic distances 
(distance modules) monotonically increases by redshift, hence any two cosmological 
models can become indistinguishable (up to a very high precision of the data) if 
we multiply the distance modules of one of them by a function form with a certain 
degree and for some particular coefficients. The Bayesian interpretation of Crossing 
Statistic is in fact hidden in these two prior assumptions that 1) /i(z) increases by 
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redshift monotonically for all cosmological models hence there are no high frequency 
fluctuations in fi(z) and 2) since the distance modules of all cosmological models 
increases by redshift, fj,(z) of any two cosmological models can become so close to 
each other at all redshifts up to an indistinguishable level if we multiply /i(z) of one 
of them to a suitable smooth function of degree n with some particular values for the 
coefficients. Coefficients of these functions are in fact the Crossing hyper-parameters 
that if all of them be simultaneously consistent to zero value then we consider the 
assumed model to be consistent to the data. 

Since the behavior of matter density and its effect on the expansion history of 
the universe is known, for any assumed model of dark energy we keep f2om as a free 
parameter along with Crossing hyper-parameters of the function form. We multiply 
the luminosity distance given by our assumed dark energy model fiDE^om, z) to the 
function Tjv(Ci, ...Cn,z) and fit this new function to the data using \ 2 minimization 
and find the best fit point in the hyper-parameter space and derive the confidence 
limits. 

Crossing function is a parametrization proposed to pick up the possible devi- 
ations from the assumed theoretical model using the data. The basic difference 
between the Bayesian interpretation of the Crossing statistic and usual parametriza- 
tions of cosmological quantities is in the fact that we are not comparing cosmological 
models within a certain framework of a parametric form and finding the best fit point 
in the parameter space and defining the confidence limits, but instead we are com- 
paring an assumed cosmological model with its own possible variations. If the data 
suggest that the deviations from the actual model must be significant we can rule 
out the model right away. In fact parameterizations of cosmological quantities can 
never be used to falsify the assumed parametric form and we can only find the best 
fit point in the parameter space while the Bayesian interpretation of Crossing statis- 
tic is designed to falsify any assumed theoretical model. Crossing function can have 
different forms but considering the shape of the data and our expectations from the 
theoretical models (once again, /i(z) should increase smoothly and monotonically by 
redshift) we expect a simple Chebyshev polynomials [22] can perform satisfactorily. 
Chebyshev polynomials have properties of orthogonality and convergence within the 
limited range of —1 < x < 1 which is appropriate in our case of study. One can 
use Gram-Schmidt process to generate a set of orthogonal basis tailored appropri- 
ately for the case of study but in our analysis and considering the range of the data 
Chebyshev polynomials are good enough. 

To falsify a theoretical model one needs to fit the given theoretical model along 
with the Crossing functions of different degrees. Looking at the data having more 
than two crossing between the actual model and the assumed model seems to be 
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unrealistic [1] and even if it is possible for some special models, it would not be de- 
tectable by the data since there are uncertainties in the data and if two models be so 
close to each other, data cannot discriminate between them. However, one needs not 
to be worry about this since assuming a Crossing function of high degree increases 
the allowed A% 2 with respect to the x 2 of the best fit point for different confidence 
levels so automatically analysis will show that by increasing the degrees of freedom 
of the function form the data would not be sensitive to pick up the tiny differences. 

In this analysis I assume Chebyshev polynomials of orders one and two as the 
Crossing functions and to show the performance of the method, I do test the spatially 
flat ACDM model for three sets of simulated data according to three dark energy 
models (all spatially flat universes), an evolving dark energy model or Kink model 
with the same parameters used in [23] 1 with Q m = 0.27, a model of dark energy 
with constant equation of state of w(z) = 0.9 with Q m = 0.27 and finally ACDM 
model with w(z) = —1 and tt m = 0.27 (correct model). Simulated data is based 
on future JDEM data with 2298 data points in the range of 0.015 < z < 1.7 [24]. 
The Kink model I used in this analysis has a special form of the equation of state of 
dark energy that at low redshifts it converges to w(z) — — 1 and at higher redshifts 
it smoothly changes to w — —0.5. This evolving equation of state of dark energy is 
not within the possibilities of CPL [25] or many other dark energy parametrizations. 
Hence using these parameterizations results to wrong reconstruction of dark energy 
if the data is based on the Kink model (look at upper-right panel of figure 6 in [23]). 

Crossing functions are given by: 

T I (C 1 ,z) = l + C 1 (— ) (2.2) 

Zmax 

and 

T n {C u C 2 , z ) = i + C 1 (—) + C 2 [2{—f - 1], (2.3) 

Zmax Zmax 

where z max is the maximum redshift in the data set. The functions that we 
fit to the data are fJ^ E {z) = fi DE (Q 0m ,z) x T^C^z) and /J% E (z) = fi DE (Q 0m ,z) x 
T n (Ci, C 2 , z) and variables are flom, Ci in case of one crossing and O 0m , Ci, C 2 in case 
of two crossings. fiDE^om, z) represents the luminosity distance of the assumed dark 

1 equation of state of dark energy in this particular model is given by: 



w(z) =w + (w m - w ) 

X 1 



l + exp(A t - 1 (l + z t )- 1 ) 
1 - cxptAr 1 ) 
exp(A t - 1 ) + exp(A t - 1 (l + zt)- 1 ) 



(2.1) 



exp(Ar 1 (l + z)- 1 ) + cxp(A t - 1 (l + zt)- 1 ). 
with the constants having the values wq = —1.0, w m = —0.5, z t = 0.5, A t = 0.05. 
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Figure 1: la and 2a confidence limits of fitting fX^j E (z) assuming the theoretical model 
being ACDM to three different simulated data sets. In the left panel the data is generated 
based on kink dark energy model with f^om = 0.27, the data in the middle panel is based 
on the ACDM model (w = —1) with £lo m = 0.27 and in the right panel data is based on 
quiescence dark energy model with w = —0.9 and £lo m = 0.27. The vertical line at C\ = 
represents the consistency of the assumed model and the data. Two horizontal lines at 
fiom = 0.25 and f2o m = 0.29 represents possible constraints we may be able to put on 
the matter density from other cosmological observations to help us falsifying cosmological 
models. 

energy model with matter density of fiom- As it was mentioned earlier one can go 
to higher orders of crossing but we will see that the simulated data we are using are 
not sensitive to the higher orders. The rest is the usual likelihood analysis based on 
X 2 statistics. In [1] we have explored that x 2 statistic is in fact equivalent to the last 
mode of Crossing statistic and in the current analysis in fact we use an approximation 
to transfer the information of the first and second Crossing modes via the Crossing 
function to the last mode which is x 2 ■ One may argue that the form of Crossing 
function may not be appropriate (to map the wrong model to the actual model) 
but considering the form and range of the data Chebyshev polynomials seems to be 
the most appropriate choice to do this mapping to the level of indistinguishability 
between the models. 

In figure. 1 we see the results of fitting jj^ E (z) assuming the theoretical model 
being ACDM to three different simulated data sets. The data in the left panel is 
based on kink dark energy model with Qo m = 0.27, the data in the middle panel is 
generated based on ACDM model (w — — 1) with f2o m = 0.27 and the right panel is 
for the data based on quiescence dark energy model with w = —0.9 and f2 0m = 0.27. 
All models are spatially flat and the data is according to the future JDEM data. 

We can Clearly see that w = —0.9 model can be ruled out with high confidence 
and the kink dark energy model shows inconsistency with the data at one to two 
sigma with no knowledge of matter density. Having a knowledge of matter density 
from other cosmological observations we can rule out the kink dark energy model 
too as it suggest the value of matter density should be about 0.32 for zero value of 
Ci for standard flat ACDM model. 
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Figure 2: la and 2a confidence limits of Ci-^om (upper panels) and C\-Ci (lower panels) 
of fitting [i'de(z) assuming the theoretical model being ACDM to three different simulated 
data sets similar to figure 1 (data is based on kink dark energy model in the left, ACDM 
in the middle and w = —0.9 in the right, all with f!o m = 0.27). The blue dots in the lower 
panels represents C\ = C2 = which shows the consistency of the assumed model and the 
data. One can see that the data are not sensitive to the second order function form Tjj to 
distinguish between the cosmological models. 



Assumed Model: ACDM 







X 2 {^de) 


x 2 {v 


11 ) 

DEI 




Data based on ACDM Model with Q 0ra = 


= 0.27 


2368.9 


236< 


3.2 


2367.9 


Data based on Kink Model with Q 0m = 


0.27 


2372.5 


236< 


3.7 


2368.5 


Data based on w = —0.9 Model with f2 0m 


= 0.27 


2375.5 


236< 


3.9 


2367.9 



Table 1: Derived \ 2 using fiDE, P^de ano - Mdb fitting the three different simulated data 
sets. Note the changes in the \ 2 by increasing the degrees of freedom of the functions. Its 
clear that in all three cases, X 2 (^de) anc ^ X 2 (a*Sb) are so c ^ ose to each other that we can 
conclude the data is not sensitive to the second order of functions. 



In figure 2 we see the results of fitting fij^(z) assuming the theoretical model 
being ACDM to three different simulated data sets similar to figure 1. 

We see that by increasing the degrees of freedom of the Crossing function, ACDM 
model becoming more consistent to all three sets of the data which shows that the 
data is not sensitive to the crossing function of second order for these particular 
models of dark energy. I should emphasize that all of these three assumed dark 
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energy models are very close to each other when we allow matter density to vary [26] . 
For instance ACDM model with Q 0m = 0.27 is extremely close to kink dark energy 
model with Qo m = 0.314. Looking at figure 2 also make it clear that there is no need 
to go to the higher orders of crossing. It basically means that the data cannot pick 
up two or three crossings of the comparing models that is clearly understandable and 
expected considering the quality and the range of the data. 

To test the sensitivity of the method to the unknown dispersion of the data, a int , 
we re-do the analysis by increasing and decreasing the intrinsic dispersion of the data 
ai n t by 20% for all data points. In figure 3 we see the results for fitting fJ^j E {z) based 
on spatially flat ACDM model to the data with under and over estimated sizes of 
the error-bars. As we see, the method still works robustly and can differ between 
the models accurately even when we have large uncertainties in estimation of the 
intrinsic dispersion of the data. 

w = —0.9 model can be ruled out with high confidence and the kink model still 
shows inconsistency with the data at one to two sigma and having a knowledge of 
matter density from other cosmological observations we can rule out the kink model 
with high confidence. 

3. Conclusion 

In summary, we have presented a new approach based on Crossing Statistics that can 
be used to distinguish between different cosmological models using Crossing functions 
and hyper-parameters. We propose that instead of parametrizing the cosmological 
quantities such as luminosity distance, Hubble parameter or equation of state of 
dark energy, we can parametrize the possible deviations from any assumed model 
and compare a model with its own variations. Performance of the method is shown 
being very much promising in distinguishing cosmological models with and without 
knowledge of matter density in comparison with other parametric and non parametric 
methods. Using this approach one can falsify the standard model (or any particular 
dark energy model) given the data without knowledge of the underlying actual model 
of the universe and its parameters. 

One may notice an interesting analogy between the first and second orders of 
Crossing functions in testing cosmological models using the cosmic distances and the 
spectral index and running of the spectral index of the primordial spectrum in CMB 
analysis. Similar to CMB analysis that we use the spectral index n s to judge about 
the consistency of the scale invariant power spectrum (Harrison-Zeldovich form of 
the primordial spectrum) to the data by looking at the derived n s and its error-bars 
(whether n s = 1 is in agreement to the CMB data), we use C\ to test the consistency 
of the assumed dark energy model to the data (whether C\ — is in agreement 
to the distance data). Similarly there is an analogy between the running of the 
spectral index in CMB studies and C*2 in our analysis. One may argue that n s is a 
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Figure 3: la and 2a confidence limits of fitting fijj E (z) assuming the theoretical model 
being ACDM to three different simulated data sets similar to figure 1 where in upper panels 
we have reduced the size of the error-bars for all data points by 20% and in the lower panels 
we have increased the size of the error-bars for all data points by 20%. We can see that 
the method is robust against the uncertainties in the data up to a very high level. 

physical parameter while C\ is not. We should note that independent of considering 
inflationary scenarios (or any physical mechanism resulting to a power-law form of 
the primordial spectrum) one can still use n s from a phenomenological perspective 
to rule out Harrison-Zeldovich form of the primordial spectrum. 

The Bayesian interpretation of the Crossing Statistic appears to us to be a 
promising method of confronting cosmological models with supernovae observations, 
and has the potential to be trivially generalized to other cosmological observations. 
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